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ABSTRACT 

A  linearized,  global  spectral  model  with  eight  levels  was  used  to  determine 
whether  the  nonlinear  interaction  between  a  planetary-scale  wave  (wavenumber 
four)  and  a  rapidly  growing  synoptic-scale  wave  (wavenumber  seven)  could 
produce  a  northeasterly  wind,  characteristic  of  East  Asian  cold  surges.  The 
amplitude  of  the  synoptic-scale  wave,  or  generic  cyclone,  was  produced  by  a 
nonlinear  Eady  model  of  the  atmosphere  that  included  friction.  The  resulting 
nonlinear  forcing  was  applied  to  either  the  first  law  of  thermodynamics,  the 
vorticity  equation,  or  both. 

The  thermal  forcing  did  not  produce  a  significant  cold  surge  response.  The 
vorticity  forcing  produced  a  respectable  cold  surge  within  48  hours.  The  results  of 
this  study  indicate  the  planetary-synoptic  wave  interaction  is  a  possible  method  for 
initiating  East  Asian  cold  surges. 
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I.    INTRODUCTION 

The  monsoon  is  a  three-dimensional,  planetary  scale  wind  regime  that  exhibits 
a  strong  seasonal  dependence.  According  to  the  Glossary  of  Meteorology 
(Huschke,  1959),  "the  primary  cause  (of  the  monsoon)  is  the  much  greater  annual 
variation  of  temperature  over  large  land  areas  compared  with  neighboring  ocean 
surfaces,  causing  an  excess  of  pressure  over  the  continents  in  the  winter  and  a 
deficit  in  summer."  The  variation  in  temperature  results  from  the  position  of  the 
sun  during  each  season.  The  shapes  of  the  continents  and  their  variable 
topographies  produce  considerable  regional  and  temporal  variability  of  monsoons. 

The  northeasterly  monsoon  that  occurs  in  East  Asia  during  the  winter  is  one  of 
the  most  energetic  circulations  of  the  atmosphere.  Even  though  the  regional 
characteristics  of  the  winter  monsoon  occur  in  East  and  Southeast  Asia,  the 
influence  on  other  components  of  circulation  can  reach  global  scales  (Lau  and 
Chang,  1987).  Because  of  the  large  scale  effects  on  the  atmosphere,  the  East  Asian 
winter  monsoon  has  been  an  area  of  active  research.  Recent  areas  of  investigation 
include  the  role  of  the  East  Asian  winter  monsoon  in  midlatitude-tropical  and  inter- 
hemispheric  interactions,  monsoonal  variations,  and  the  forcing  mechanisms 
responsible  for  small  scale  monsoonal  variations. 

The  East  Asian  winter  monsoon  is  associated  with  the  thermally  direct  Hadley 
circulation,  or  cell,  that  occurs  over  the  area  in  winter.  The  ascending  branch  of 
the  Hadley  cell,  and  major  monsoonal  convective  zone,  migrates  from  its 
summertime  position  over  India  to  the  maritime  continent  of  Borneo/Indonesia 
(Ramage,  1971).  The  latent  heat  released  in  the  upper  levels,  due  to  intense 
convection  over  the  maritime  continent,  is  transferred  poleward.     The  cold 


Siberian  high  pressure  system,  or  anticyclone,  is  the  heat  sink  area  for  the  upper- 
level  poleward  moving  warm  air.  The  Siberian  anticyclone,  in  conjunction  with 
the  descending  branch  of  the  Hadley  cell,  produces  a  large  area  of  subsidence,  and 
thus  dominates  the  Southeast  Asian  winter  (Ramage,  1971).  The  Himalayas  block 
the  southward  movement  of  extremely  cold  surface  air  from  the  Siberian 
anticyclone.  The  only  effective  outflow  region  is  to  the  southeast. 

Boyle  and  Chen  (1987)  documented  the  wintertime  surface  and  20  kPa  wind 
fields  for  the  period  1973  to  1984.  The  predominant  northeasterly  flow,  at  the 
surface,  associated  with  the  East  Asian  winter  monsoon  is  shown  in  Fig.  1.1.  The 
blacked  out  regions  correspond  to  terrain  heights  above  1000  m.  A  vector  length 
of  5°  longitude  corresponds  to  a  ten  meters  per  second  wind,  and  the  isotach 
contour  interval  is  2.5  m-s-1.  The  intense  subtropical  jet  stream  over  Southeast 
Asia,  shown  in  Fig.  1.2,  is  caused  by  the  intense  baroclinic  zone  between  the  warm 
tropics  and  frigid  Siberian  area.  A  wind  vector  length  of  5°  corresponds  to  a  100 
m-s-1  wind,  and  the  isotach  interval  is  10  ms_1.  The  climatological  averaged 
January  50  kPa  geopotential  height  and  wind  field  (Fig.  1.3)  shows  a  dominant  long 
wave  trough  centered  over  the  Sea  of  Japan..  The  contour  interval  is  80  m. 
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Figure  1.1.    Average  Wintertime  Surface  Wind  (Boyle  and  Chen, 

1987) 
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Figure  1.2.    Average  Wintertime  20  kPa  Wind  (Boyle  and  Chen,  1987) 
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Figure  1.3.    Climatological  January  50  kPa  Height  Field  (Boyle  and 

Chen,   1987) 

The  Hadley  circulation,  and  the  effect  of  the  Himalayas,  set  up  the  planetary 

circulation  of  the  winter  monsoon.  However,  Boyle  and  Chen  (1987)  indicate  that 

transient  synoptic-scale  waves  shape  the  final  form  of  the  Siberian  High.    As 

synoptic  waves  propagate  along  the  longwave  trough,  surface  cyclones  and 

anticyclones  develop  due  to  the  intense  baroclinicity.  Cyclones  usually  develop  off 

the  west  coast  of  Japan  in  the  area  of  strong  upper-level  positive  vorticity  advection 

and  low-level  warm  air  advection,  and  track  to  the  northeast.  Anticyclones  develop 

near  the  southern  extent  of  the  Siberian  High,  due  to  the  upper-level  negative 

vorticity  advection  and  low-level  cold  air  advection,  and  track  toward  the 

southwest  over  China.  When  the  pressure  gradient  between  the  China  anticyclone 

and  the  cyclone  off  Japan  tightens  rapidly,  significant  ageostrophic  motion  results. 

The  cross-isobaric  ageostrophic  flow  accelerates  toward  lower  pressure  near  the 

cyclone  to  the  east  and  to  the  Intertropical  Convergence  Zone  (ITCZ)  to  the  south. 

The  northeasterly  cold  wind  is  the  "cold  surge"  in  the  winter  monsoon.  Cold  surges 


usually  reach  the  equatorial  South  China  Sea  in  12  to  24  hours.  The  enhanced 
northerly  flow  intensifies  the  tropical  convection  over  the  maritime  continent,  and 
thus  strengthens  the  Hadley  circulation.  The  cold  surge  ends  when  the  midlatitude 
trough-ridge  pattern  moves  far  enough  east  to  diminish  the  China  anticyclone,  and 
consequently  the  pressure  gradient.  Boyle  and  Chen  (1987)  emphasize  that  cold 
surges  are  caused  by  the  interaction  of  synoptic  and  planetary  waves.  They 
conclude  that  cold  surges  are  dynamically  forced,  and  thus  must  be  considered 
separate  entities  from  the  Siberian  anticyclone. 

The  East  Asian  cold  surge  has  two  stages.  They  are  separated  by  a  few  hours  to 
one  day,  depending  on  the  location  of  the  observing  site  (Lim  and  Chang,  1981; 
Chang  et  al.,  1983).  The  first  stage  is  the  pressure  surge.  It  is  the  leading  edge  of 
the  air  accelerating  towards  the  equator.  The  pressure  surge  propagates  with  the 
speed  of  internal  gravity  waves.  The  second  stage  is  a  frontal  passage  that  moves 
with  advective  speeds.  It  is  defined  by  a  sharp  decrease  in  surface  temperature  and 
dew-point  temperature. 

The  onset  of  the  East  Asian  cold  surges  are  defined  in  many  ways  (Boyle  and 
Chen,  1987).  The  three  most  common  definitions  are:  1)  A  drop  in  surface 
temperature  at  Hong  Kong  of  five  degrees  Celsius,  or  more;  2)  An  increase  of  the 
surface  pressure  gradient  between  coastal  and  central  China  of  at  least  0.5  kPa;  and 
3)  A  prevalent  northerly  surface  flow  over  the  South  China  Sea  with  speeds 
exceeding  five  meters  per  second(Lau  and  Chang,  1987). 

The  cause  or  nature  of  East  Asian  cold  surges  has  been  examined  by  Lim  and 
Chang  (1981),  Baker  (1983),  Bashford  (1985)  and  Harris  (1985).  Lim  and 
Chang  (1981)  used  linearized  shallow- water  equations  on  an  equatorial  beta-plane 
to  simulate  the  response  of  the  tropics  to  a  midlatitude  pressure  surge.  They  did  not 


include  planetary  boundary  layer  friction  or  orography.  Lim  and  Chang  (1981) 
found  that  synoptic  scale  forcing  in  the  midlatitudes  produced  Rossby-type  waves 
that  propagated  into  the  tropics.  The  northeast-southwest  tilt  in  the  pressure  field, 
typical  of  Rossby  waves,  and  the  northeasterly  flow  similar  to  cold  surge  events  is 
shown  in  Fig.  1.4. 
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Figure  1.4.    Wind  Velocity  Vectors  for  a  Barotropic  Model  (Lim  and 

Chang,  1981) 

Baker  (1983)  used  a  global,  six-layer,  primitive  equation  model  to  examine  the 
interaction  of  a  midlatitude  baroclinic  wave  with  topography.  Baker's  results 
indicated  that  well  developed  baroclinic  waves  could  initiate  a  cold  surges,  but  the 
surges  were  limited,  and  weak.  Baker  concluded  that  other  forcing  mechanisms 
were  required  to  simulate  cold  surges. 

Bashford  (1985)  used  an  eight-layer,  spectral,  primitive  equation  model  with 
an  analytical  heat  source  to  study  the  effects  of  planetary  scale  motion  on  cold 
surges.  The  heat  source  function  of  Bashford's  baroclinic  model  were  similar  to 
that  used  by  Lim  and  Chang  (1981)  in  their  barotropic  model.  Bashford  found  that 
a  planetary  wave  (wavenumber  three)  with  a  deep  thermal  forcing  could  produce  a 
cold  surge  response. 


Harris  (1985)  used  the  same  model  as  Baker  (1983)  to  study  the  interaction  of  a 
baroclinic  wave  (wavenumber  eight)  with  a  planetary  wave  (wavenumber  four). 
Harris  found  the  synoptic  wave  alone  did  not  produce  a  cold  surge,  but  the 
synoptic-planetary-wave  interaction  produced  a  significant  cold  surge  response. 
Harris'  results  may  be  tainted  because  the  planetary  wave  extended  unrealistically 
far  south. 

This  study  will  use  a  global,  eight-layer,  primitive  equation  model,  and  will  be 
conducted  in  three  phases.  The  first  phase  will  examine  the  effect  of  a  forcing  with 
a  single  frequency  derived  from  a  forced  Rossby  wave,  midlatitude  beta-plane 
solution.  The  second  phase  will  repeat  Bashford's  (1985)  experiment  with  an 
impulse  forcing  that  includes  a  spectrum  of  frequencies.  This  effort  will  attempt  to 
reproduce  the  results  of  Lim  and  Chang  (1981).  The  third  phase  will  examine  the 
interaction  of  a  planetary  wave  (wavenumber  four)  with  a  synoptic  wave 
(wavenumber  seven).  The  synoptic  wave,  or  generic  cyclone,  will  be  produced  by 
a  nonlinear  Eady  model  of  the  atmosphere  (Peng,  1982).  The  basic  state  or  mean 
flow  of  the  planetary  wave,  and  the  phase  speed  of  the  synoptic  wave,  are  both  equal 
to  zero.  Instead  of  combining  the  two  waves  together  within  the  global  model,  as 
done  by  Harris  (1985),  this  study  will  compute  the  nonlinear  interactions 
analytically.  The  resulting  wave  (wavenumber  three)  is  used  in  the  linear,  global 
model.  To  aid  in  the  analysis  of  each  experiment,  the  solutions  from  the  global 
model  will  be  projected  onto  the  vertical  normal  modes. 


II.     SIMPLE  BETA-PLANE  SOLUTION 

Since  the  baroclinic,  primitive  equation  model  used  in  this  study  is  linear,  and 
the  mean  or  basic  flow  is  zero,  the  full  baroclinic  atmosphere  can  be  represented  by 
the  sum  of  vertical  normal  modes.  The  behavior  of  each  vertical  mode  solution  is 
similar  to  a  barotropic  shallow  water  system  with  the  appropriate  mean,  or 
equivalent,  depth  (Lim  and  Chang,  1987).  In  other  words,  each  mode  of  the 
baroclinic  system  will  behave  like  a  shallow  water  system.  To  simulate  the 
variation  in  the  Coriolis  parameter,  and  simplify  the  mathematics,  the  system  of 
shallow  water  equations  is  solved  on  a  midlatitude  beta-plane.  The  validity  of  using 
a  midlatitude  beta-plane  is  discussed  by  Lindzen  (1967). 

The  scaled  shallow  water  equations  on  the  midlatitude  beta-plane,  with  a 
forcing  function  added  to  the  continuity  equation,  are: 

dt  fo  (2.1) 

^+vrv(£  +  poy)+fbD=0 

ol  (2.2) 

vV-foC=o  (23) 

where: 

<j)'         -  perturbation  in  the  geopotential  height  field 

\\f         -  geostrophic  streamfunction 

H         -  mean  (basic)  height  field  or  equivalent  depth 

D         -  divergence 

V         -  two  dimensional  wind  velocity  vector 


£  -  vorticity 

fo         -  Coriolis  parameter  (constant) 

t  -  time 

Po        -  gradient  of  Coriolis  parameter  in  the  north-south  direction 

y  -  north-south  position 

S  -  forcing  function 

If  the  vorticity  is  assumed  geostrophic,  then  £  =  V2\j/,  and  Eq.  2.3  yields  \|/  = 
4>'/foor  0'  =  fo\|/.  When  <j)'  is  introduced  into  Eqs.  2.1  and  2.2,  and  the  nonlinear 
terms  are  eliminated,  Eqs.  2.1  and  2.2  become: 

f(^  +  gHD=-Mis 

&  f0  (2.4) 

^  +  po|^  +  foD=0 

<*         <**  (2.5) 

When  the  divergence  is  eliminated  between  Eqs.  2.4  and  2.5,  the 
quasigeostrophic  potential  vorticity  equation  results: 

dt\      gH    /  dx  (2.6) 

Eq.  2.6  is  solved  by  writing  the  streamfunction  and  forcing  function  as  follows: 
¥=vF(y)e'-(Kx-cDt)  (2?) 

S=/s(y)eJ'(Kx-wO  (2.8) 

where  K  is  the  dimensional  wavenumber,  and  (0  is  the  frequency.  The  latitudinal 
structure  of  the  forcing  function  is: 


[Wgr]    lyl<w 
s(y)= 

0  |y|>W 

1  y  (2.9) 

where  So  is  the  magnitude  of  the  forcing,  and  W  is  half  of  the  width  of  the  forcing 
in  the  north-south  direction. 

When  Eqs.  2.7  and  2.8  are  substituted  into  Eq.  2.6,  the  following    ordinary 
differential  equation  results: 

d2vr  ia2y^  s(y) 

dy2  w  (2.10) 

where: 

^     w  (2.11) 

or: 

-Kp0 


00  = 


K2  +  a2  +  -S- 

^  (2.12) 

The  homogeneous  solution  to  Eq.  2.10  is: 

4>(y)  =  e±<ay  (2.13) 

Since  the  latitudinal  structure  is  symmetrical  about  y  =  0,  Eq.  2.10  will  only  be 
solved  for  lyl  <  -W  and  y  <  -W.  When  y  <  -W,  the  negative  exponent  of  Eq.  2.13 
is  used,  and  when  y  >  W  the  positive  exponent  is  used.  If  a2  >  0,  the  particular 
solution  to  Eq  2. 1 0  for  the  interval  lyl  <  W  is: 


i}/(v)=_So_  +  So vw 

2a2     2a2/_7Lj2 


W/  (2.14) 
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When  Eqs.  2.13  and  2.14  are  combined,  the  general  solution  (with  a2  >  0)  is: 


¥(y)=Sl 


1 


COS 


Iwl 


a2-(^f 
W/ 


xF(y)  =  C2e-'ay 


y<-W 


+  Ci(e''ay  +  e-''oy)         iyi<w 

(2.15) 

(2.16) 

Since  a2  >  0,  a  can  be  positive  or  negative.  When  a  >  0,  the  solution  to  Eq.  2.16 
will  have  a  north-south  phase  structure  that  tilts  from  the  northeast  to  the  southwest 
for  y  <  -W.  A  northeast-southwest  tilt  in  phase  structure  will  propagate  energy 
away  from  the  source  region.  This  can  be  shown  by  applying  a  radiation  condition 
which  requires  that  the  wave  have  a  group  velocity  moving  away  from  the  source. 
If  a  <  0,  the  solution  to  Eq.  2.16  will  have  a  phase  structure  that  tilts  from  the 
northwest  to  the  southeast,  and  will  propagate  energy  toward  the  source  region. 
For  the  region  y  >  W,  a  <  0  will  produce  a  northwest-southeast  tilt  in  the  phase 
structure,  and  propagate  energy  northward  away  from  the  source  region. 

The  two  constants  are  determined  by  equating  Eqs.  2.15  and  2.16,  and  also 
equating  their  derivatives,  at  y  =  -W.  When  Eqs.  2.15  and  2.16  are  substituted  into 
the  real  part  of  Eq.  2.10,  the  complete  propagating  solution  is: 


/ 


xj/_So 

"  2 


1 


cos 


+  ■ 


Tcy\  \ 

wi 


a2.fjL 
W 


1 


cos(Kx-cot)- 


a2-W 
W' 


¥  =  ^sin(aW)/-J 1 


or 


or 


;  (ay)oos  (Kx  +  aW  -  cot) 


sin(Kx  +  ay-cot) 


W'  / 


lyl<W 


y<-W 


(2.17) 


(2.18) 
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If  a2  <  0  and  a2  =  -{I2,  Eq.  2.10  becomes: 

d2vF    u2^=   s(y) 
dy2  w 

and  Eq.  2.1 1  becomes: 

gl        CO 

The  homogeneous  solution  of  Eq.  2.19  is: 
and  the  particular  solution  for  lyl  <  W  is: 

VJ/(y)  =  _J>0__S0_  \W/ 


(2.19) 


(2.20) 


(2.21) 


4*    V+(^f 


(2.22) 


When  y  <  -W,  the  positive  exponent  of  Eq.  2.21  is  used. 

When  Eqs.  2.21  and  2.22  are  combined,  the  general  solution  (with  a2  <  0  or 
H2  >  0)  is: 


<F(y)  = 


=-Sq 


cos  - 


/7CV\    \ 


1  +       \w/ 

^2     n?+(JEf 

\       K     \w/ 


w/  / 


+  Ci(eW  +  e-w)  lyl<W 

(2.23) 

^(y)  =  C2eW  y<_W  (2.24) 

Since  |l2  >  0,  |i  can  be  positive  or  negative.  When  \x  >  0,  the  solution  to  Eq.  2.24 
will  trap  energy  near  the  source.  If  |i<  0,  then  the  solution  would  grow  to  infinity 
for  large  lyl. 
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The  constants  can  be  solved  in  the  same  manner  as  used  to  get  Eqs.  2.17  and 
2.18.  When  Eqs.  2.23  and  2.24  are  substituted  into  the  real  part  of  Eq.  2.10,  the 
complete  trapped  solution  is: 


4>=So  I  __1 


cos  — 

lw/ 


2[_»+(JLf4 


lw/ 


e-^wcosh(jiy; 


*  Htf 


i] 


cos  (Kx  -  cot) 


¥=^sinh(^Wp 1 Wysin(Kx-cot) 


lw/  / 


lyl<  W 


y<-W 


(2.25) 


(2.26) 

The  derivation  above  indicates  that  under  certain  conditions,  or  frequencies, 
the  solution  to  the  shallow  water  equations  can  either  propagate  energy  away  from 
the  source  region,  or  trap  energy  near  the  source.  The  experiments  in  Phase  I  will 
examine  the  effects  of  forcing  the  primitive  equation  model  with  a  variety  of 
frequencies. 
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III.    MODEL  DESCRIPTION 

The  model  used  in  this  study  is  a  baroclinic  spectral  transform  primitive 
equation  model  as  described  by  Haltiner  and  Williams  (1980),  and  used  by  McAtee 
(1984).  The  model  discussion  of  this  chapter  is  taken  directly  from  McAtee.  The 
model  is  configured  to  include  friction,  diabatic  heating  and  a  vorticity  forcing. 
The  specifics  of  how  diabatic  heating  and  vorticity  are  included  in  the  model  will 
be  discussed  in  the  next  section.  Friction  is  not  included  in  this  study.  The  basic 
equations  of  the  model,  in  sigma  coordinates,  are  as  follows: 

^  =  -V-(;+f)v-k-Vx(RWq  +  o^|  +  k-VxF+A 
at  \  do  J 


—  =  kVx(C  +  f)V-V-(RTv7q  +  o— I- v1<{>+^-)  + V-F 

at  ^    '      \     4     del     v    2 1 


^  =  -D-VVq-3a 


(3.1) 

(3.2) 
&  do  (3.3) 

(3.4) 

(3.5) 


^=-vve- 

0—  +  Q 

do 

^  _    RT 
do       o 

where: 

c     - 

vorticity 

D 

divergence 

T 

temperature 

t 

time 
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0  -  potential  temperature 

n  -  surface  pressure 

V  -  horizontal  velocity  vector 

<J>  -  geopotential  height 

R  -  gas  constant 

Cp  -  specific  heat  at  constant  pressure 

f  -  Coriolis  parameter 

P  -  PT 

o         -  vertical  coordinate  (a  = ) 

n-pT 

da 

a         -  vertical  velocity  (6  = -r-) 

q         -  ln(P) 

P         -  pressure 

Pt        -  pressure  at  the  top  of  the  model 

F         -  frictional  force 

Q         -  diabatic  heat  forcing  function 

A         -  vorticity  forcing  function 

The  continuity  equation  (3.3)  is  rewritten  by  integrating  with  respect  to  sigma. 
and  applying  the  boundary  conditions  a  (0)  =  a  (1)  =  0.0.  Thus  the  integral  of 
Eq.  3.3  may  be  written: 

&  (3.6) 


where: 
(")=  j   (  )do 

Jo 
G  =  V-Vq 
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The  vertical  velocity,  6,  is  derived  diagnostically  by  substituting  Eq.  3.6  into 
Eq.  3.3,  integrating  in  the  vertical,  and  using  a  (0)  =  0. 

o=(D  +  G)a-      (D  +  G)do 

Jo  (3.7) 

The  first  law  of  thermodynamics,  Eq.  3.4,  is  written: 

?I  =  -VVT-aKd^!^  +  KT(^.  +  VVq]+Q 
9t  Be  \8t  }  (38) 

where  k  equals  R/Cp.  To  apply  semi-implicit  differencing,  it  is  necessary  to  divide 

the  temperature  into  parts  as  follows: 

T  =  T*(o)  +  T(oXw)  (3.9) 

where: 

T*        -  appropriately  averaged  temperature 

T'        -  perturbation  temperature 

X         -  longitude 

cp         -  latitude 

The  basic  equations  are  conveniently  written  in  spherical  coordinates  by 
defining  the  following  operator: 

cos2(p3^    coscpatp  (31Q) 

Using  Eqs.  3.9  and  3.10,  the  basic  equations  can  be  written  as  follows: 

3C 


■  =  -<x(A,B)+A 


oD 


=  -a(A3)-VlE  +  <j)+RT*q) 
at 


(3.11) 
(3.12) 


16 


—  =  -a(uT,vr)+DT'-oKd3^^  +  KT(G-G-D)  +  Q 
*  d°  (3.13) 


|=-(D+G) 


(3.14) 


a^-=-RT 

3o  (3.15) 

where: 

a     ir    At     -dV    RT      7    *J   coscp 

A  ^  £ + f  )U + o— + -— cosfy -^- -  — — F<p 

3o     r2  ax,      r      Y 

B=(C+f)v-a!L^+^ 

ao  r2  ax,     r 
G:r   u   dq ,  v  i    dq 

cos29  3^       coscp  3(p 

E=u2  +  v2 


u= 
v= 


2cos-cp 

ucoscp 
r 

vcoscp 


and: 

u  -  zonal  component  of  the  velocity  vector  (V) 

v  -  meridional  component  of  the  velocity  vector  (V) 

r  -  mean  radius  of  the  Earth 

Eqs.  3.10  through  3.15  are  the  basic  equations  used  in  the  model.  These  equations 

are  represented  spectrally  in  the  horizontal  and  finite  differenced  in  the  vertical. 
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A.   VERTICAL  STRUCTURE 

The  vertical  structure  of  the  model  follows  the  development  given  by  Arakawa 
and  Suarez  (1983).  The  variables  are  staggered  in  a  so  that  £,  D,  U,  V  and  T  are 
carried  at  the  mid-point  of  each  layer,  where  a  =dk.  The  variable  a  is  carried  at  the 
top  and  bottom  of  each  layer,  where  o  =  Ok-  The  vertical  structure  is  illustrated  in 
Fig.  3.1.  The  finite  difference  form  of  Eqs.  3.10  through  3.15  are: 


-^  =  -ak(A3)  +  Ak 
at 

^  =  ak(B,-A)-Vfy>k  +  RT*q  +  Ek) 
at 


Pk         Pk+ 

4>LM  =  <))S  +  CpTLN1 


(3.16) 
(3.17) 
(3.18) 


at       l  ' 

^=-ak(ur,vr)r  +  T;(Gk-G-D)- 

at 

-Uak+iBii-^k-Tk+i  -  Tk)  +  akAk.!(Tk  -  -Sl-ThIUq 

A(V  UTc+1  I  »         Pk-1         I)  (3.19) 

0k-<l>k+l  =  C^Pk+i-Pk)(^+^Ik±I 


(3.20) 


Plm     /  (3.21) 

k 

6k+i=ok+i(G+D)-X(Gj  +  Dj)Aaj 

H  (3.22) 

where: 
^(cfcuIlF 

p=U±KJ 

Ok+1  -  C?k 
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=(c  i  flu  i  °k^Vk+1  -  vk)+<7k(vk  -  yk-i)  t  (gjjggM    i    fo  (oosylp 

|       r2        /coscpacp   I    r    F<P 


2Aak 


B  =  (C+f)v   ta^i(Uk+i-Uk)+gk(Uk-Uk.i)   (CpRTk|    }     3g  t(cos<pL 
lS     '  JAa,  {      r      'coscp  do    l    r    ^ 


G  = 


2Aak  r      '«*<P3cp 

vk-vn 


n 

D  =  VVk 

Ak=  ^-Pk 
Pk+i-Pk 

B      Pk-n-Pk  =  1.A 
Pk+i-Pk 
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P  =  20.0kPa  _     o,   a,   0,      o  =  0.000 

—  —    — Ci  D,  T,  Ac,  0, 

AAA  #        A 

—  C2  D2  T2  °2  °2  02 0  =  0.125 

£2  ^2  ~^2  AO2  02 

-C3  D3  T3  a3  a3  J3  a  =  0.250 

—  —     —     —     —       C3  D3  T3  A<*3  03 

-£4  ^  T4  a4  a4  04 a  =  0.375 

_____     C^  D4  T^  Aa*  0^ 

-C5  D5  T5  a5  a5  05  a  =  0.500 

—     —     —       fc  D5  T5  Aa5  05 

-Ce  D6  T6  a6  o5  0e  a  =  0.625 

—  —     —     —     —       Ce   D6  T6  A^6  06 

-£7  D?  T7  o7  a7  07  °  =  °-750 

—  —    —    —    —     C?  D7  T7  Aa7  07 

-Ce  D8  T8  a8  o8  08  °  =  0.875 

—  —   —   —   —     Ce  D6  Te  Ao8  08 

P5  ■  1 01 .325  kPa  a  ,       77777m^  a  ■  1 .000 


Figure  3.1.     Vertical  Structure  of  the  Primitive  Equation  Model 

Eqs.  3.16  through  3.22  are  written  in  matrix  form,  so  the  terms  on  the  right 
hand  sides  contain  all  the  terms  that  are  to  be  evaluated  explicitly,  and  the  terms  on 
the  left  hand  sides  contain  those  terms  that  are  to  be  evaluated  implicitly.  Eq.  3.20 
is  combined  with  an  integrated  finite  difference  form  of  Eq.  3.15  to  obtain: 

<J>=|C|T  +  <|)S  (3.23) 

where  ICI  is  a  square  matrix  and  the  other  quantities  are  column  vectors.  The  finite 
difference  form  of  the  surface  pressure  tendency  equation,  Eq.  3.14,  is  : 

<*      k  (3.24) 

which  can  be  written  in  matrix  form  as: 
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^.  =  -NT(G  +  D) 

«  (3.25) 

where  NT  is  the  transpose  of  a  constant  column  vector.    Similarly,  Eq.  3.22  is 

written: 

o=Z(G  +  D)  (3.26) 

The  next  to  last  term  in  Eq.  3.19  is: 
^k+iBj— -Mk+i  -  Tk+okAk.]Tk  -  — M^i 


\Pk+]  /  \       Pk-i        /  (3.27) 

For  the  purpose  of  semi-implicit  formulation,  the  temperature  is  separated 
according  to  Eq.  3.9.  The  mean  part  of  that  term  is  written: 


iPk+1 


+dkAk-iTk-^-Tk.1=|M|(G+D) 

I       Pk-i        /  (3.28) 


Eqs.  3.17,  3.19,  3.23  and  3.25  are  now  written: 

—  +  V2(({)    /  +  |R|T*q)=a(B,-A)-V2E  =  KD 
dt 

^3-  +  NTD  =  -NTG 
dt 


^+|q|d  =  kt 

dt 


(3.29) 
(3.30) 


(3.31) 

<K=|c|t  (332) 

where  IQI  =  IMI  +  kI^N7"  and  $'  =  <})  -  <J)s-  The  variables  Kd  and  Kt  represent  terms 
that  are  not  explicitly  separated  out. 

The  semi-implicit  time  differencing  is  achieved  by  evaluating  the  terms  on  the 
left  hand  sides  of  Eqs.  3.29  through  3.31  implicitly.    The  remaining  terms  and 
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Eq.  3.32  are  evaluated  explicitly  using  leapfrog  differencing.    The  difference 

equations  are  written: 

D„+i+AtV2(crrn+i+|R|T*qn+1)= 

Dn.1-AtV2flC|Tn-i+|R|T*qn.1)  +  2At(KD)n  a33) 

Tn+i+At|Q|Dn+i=Tn.i-At|Q|Dn.i-2At(KT)n  (3.34) 

qn+1  +  AtNTD  =  qI1.1-AtNTDn.1-2AtNTGn  (3.35) 

Now,  the  following  equation  for  D  is  found  by  substituting  Eqs.  3.34  and  3.35  into 
Eq.  3.33: 

BJDn+1  =  B+Dn-i  +  2At(KD)n  -  2AtV^  C|  Tn.j  + 1  R|  1\.,  + 1  C|  (KT)n)- 

|R|T*NTGn  (3.36) 

where  the  matrix  operator  B  is: 

B±=At^|C||Q|±|R|TNT)V2  +  |l|  (3>37) 

and  III  is  the  identity  matrix. 

B.    SPECTRAL  FORMULATION 

The  equations  3.10  through  3.15  are  represented  spectrally  in  the  horizontal. 
The  variables  are  represented  as  follows: 

c(M>,C,t)=S    t   CH1(a,«(sin(p)e^=i    £    C^Y? 

m=-J  n=Hm|  m=-J  «Hm|  (3.38) 

where  C  is  some  variable,  and: 

c™(H-ircrim 
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where  m  is  the  zonal  wavenumber,  n  is  the  meridional  index,  and  n  -  Iml  gives  the 
number  of  zeros  between  the  poles  (-1  <  sin  9  <  1)  of  the  associated  Legendre 
function.  Triangular  truncation  is  used  in  this  study,  with  the  truncation  limit,  J, 
equal  to  three.  The  non  dimensional  zonal  coordinate  index,  X,  equals  (s  -  l)/2 
where  1  <  s  <  16.  Note  that  the  separation  is  such  that  the  coefficients  C™  are 
functions  of  time  and  the  vertical  and  spherical  harmonic  Y£  are  horizontal 
functions  of  space.  The  normalization  and  orthogonal  properties  of  Y„  allow  the 
coefficients  to  be  obtained  as  follows: 

qr=(C,YlT}^^-|       J     CY™dcpdX 

KJo    J*  (3.39) 

The  nonlinear  terms  are  computed  using  the  transform  method  following  Haltiner 
and  Williams  (1980).  The  longitudinal  direction  is  treated  with  a  Fast  Fourier 
Transform  and  the  latitudinal  direction  uses  Gaussian  Quadrature.  The  number  of 
latitudes,  N,  and  longitudes,  M,  satisfy:  N  >  3J/2  +  1  and  M  >  3(J  -  1)  +  1.  For  this 
study,  N  =  60  and  M  =  48. 
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IV.    FORCING  FUNCTIONS 

Chapter  III  of  this  paper  describes  the  global  spectral  model.  The  vorticity 
equation  (3.1)  and  the  thermodynamic  equation  (3.4)  use  forcing  functions,  A  and 
Q,  respectively.  This  chapter  describes  the  forcing  functions. 

A.   PHASE  I  -  ROSSBY  WAVE  FORCING 

The  thermal  and  vorticity  forcing  of  Phase  I  follows  the  development  of 

Chapter  II.  The  forcing  function  is: 

F=(Fa)(Fs)(F9)(Fr)  (41) 

where  Fa  is  the  amplitude  of  either  the  thermal  or  vorticity  forcing.   The  vertical 

structure  (Fs),  shown  in  Fig.  4.1,  is: 

(l.O-FAo)sinh2GLQ2  _>r_ 

sinhM.O 

tanhV^  0<(W 

tanhl.O  (4.2) 


Fs  = 


where  Fao  is  the  amplitude  of  either  forcing  at  the  surface  (for  this  study  Fao  =  0), 

and: 

Glq=  q-°max- 


<7HI  = 


1.0-Omax  (4.3) 


<5max  "  °min  (4.4) 
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FORCING 


Figure  4.1.    Vertical  Structure  of  Forced  Rossby  Wave  Solution 

The  meridional  structure  (F^)  is: 

F  -a^fr-^)) 

9  \Yn-Ys/  (4.5) 

where  <p  is  the  latitude,  and  Yn  -  Ys  is  the  north-south  period  of  forcing. 

The  development  in  Chapter  II  uses  a  cos2  latitudinal  function,  where  the  actual 
latitudinal  forcing  of  the  model  is  sin2.  The  difference  is  insignificant.  The  forcing 
is  equal  to  zero  in  the  Southern  Hemisphere,  and  outside  of  the  latitudinal  limits  in 
the  Northern  Hemisphere  (Fig.  4.2). 
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Figure  4.2.     Latitudinal  Structure  of  Forced  Rossby  Wave  Solution 
The  amplitude,  vertical  and  latitudinal  components  of  the  forcing  function 
discussed  so  far  are  similar  to  Bashford  (1985).  Both  the  temporal  and  longitudinal 
structure  of  the  forcing  functions  are  included  in  the  function  Fr  shown  below: 

-,          (2jt(X-Xw)   ,J 
Fr  =  cos| — ■ —  -cot  I 


Xe-Xw 


(4.6) 


where  Xw  -  Xe  is  the  East  -  West  period  of  forcing,  and  X  is  the  longitude.  The 
frequency  (co)  is  given  by  Eq.  2.11. 

B.    PHASE  II  -  IMPULSE  FORCING 

Bashford  (1985)  used  an  impulse  forcing  in  the  thermal  equation  similar  to 
Lim  and  Chang  (1981).  In  this  study  the  impulse  is  also  added  to  the  vorticity 
equation.  Since  this  impulse  forcing  includes  a  spectrum  frequencies,  it  is  more 
realistic  than  the  single-frequency  forcing  described  in  Phase  I.  The  impulse 
forcing  is  given  by: 
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F=(FA)(Fs)(F(p)(Fx)(FT) 


(4.7) 


where  Fa,  Fs,  and  F^  are  the  same  as  Phase  I.  The  longitudinal  part  (Fx)  is  the 
same  as  Phase  I,  but  it  is  separate  from  the  temporal  function.  The  impulse  is 
generated  by  the  temporal  function  (Fj): 


_  t2   -1 


Fr=-Crex 

2x3 


(4.8) 


where  t  is  time  and  x  is  the  time  scale.  The  impulse  peaks  at  t  =  2x  and  then  decays 
exponentially  to  zero.  The  solutions  in  this  study  use  i  =  .5  day,  so  the  peak  of  the 
impulse  occurs  at  one  day  (Fig.4.3). 


Figure  4.3.    Temporal  Structure  of  the  Impulse  Forcing 
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C.    PHASE  III  -  GENERIC  CYCLONE  FORCING 
1.    Thermal   Forcing 

The  thermal  forcing  from  the  generic  cyclone  represents  the  nonlinear 
advection  of  the  synoptic  temperature  field  by  planetary  waves,  and  vice  versa. 
The  nonlinear  advection  terms  are  given  below: 
Q  =  -  Vl-VTs  -  Vs-VTl  (49) 

where: 

Vl       -  two  dimensional  planetary  wind  vector 

Vs       -  two  dimensional  synoptic  wind  vector 

Tl       -  planetary  scale  temperature  field 

Ts        -  synoptic  scale  temperature  field 

Eq.  4.9  is  written  in  scalar  form  as  follows: 

Q.-Ul^-Vl^-Us^-Vs^ 

dx  dy  dx  dy  (410) 

The  quasistationary  geopotential  height  field  for  the  planetary  wave  is: 

<j)L  =  A(P,t)W(y)cos(KLx)  (4]1) 

where  A(P,t)  is  the  amplitude  factor  for  the  planetary  wave  (constant  for  this 
study),  and  Kl  is  the  planetary  wave  number  (for  this  study  Kl  =  4).  The  latitudinal 
structure  of  the  planetary  wave  is: 

W(y)  =  sin(^^l| 

UPN-Yps/  (4.12) 

where  Yps  -  Ypn  is  the  north-south  period  of  forcing  for  the  planetary  wave. 

The  geopotential  height  field  for  the  synoptic  wave  is: 
0s  =  N(y)(a(P,t)oos(Ksx-vt)+b(P>t)sin(Ksx-vt))  (4  13) 
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where  a(P,t)  and  b(P,t)  are  the  vertical  amplitudes  of  the  synoptic  scale  wave.  The 
functions,  a(P,t)  and  b(P,t),  represent  a  baroclinic  wave  generated  from  an  Eady 
model,  to  be  discussed  later.  The  value  K$  is  the  synoptic  wavenumber  (for  this 
study  Ks  =  7),  and  v  is  the  frequency.  The  latitudinal  structure  of  the  synoptic-scale 
wave  is: 

My-YssH 


N(y)=sin 


\YsN-Yss/  (4.14) 

where  Ysn  -  Yss  is  the  north-south  period  of  forcing  for  the  synoptic  wave. 
The  ideal  gas  law: 

P  =  PRT  (4.15) 


and  the  hydrostatic  equation: 

aT'pg 


are  combined  to  yield: 


(4.16) 


Ts=.p^=.pJ^s(Ksx.vt)+absin(Ksx.vt) 
Rap     r  ap  ap 


(4.17) 


and: 


RdP       R8P  (4.lg) 

The  planetary  wave  amplitude  function,  A(P,t),  is  constant  for  this  study 
so  Tl  =  0  and  Eq.  4.10  simplifies  to: 

Q  =  -UL^-VL^J> 

dx  ay  (4  19) 

The  planetary  and  synoptic  vorticity  fields  are: 
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k=iv>L=^ 


/   2  \ 

'aw    2    ^ 

K£W 


idy: 


cos(Klx) 


/ 


(4.20) 


Cs=4v2*s^ 


f 


f\9y2 


-KsN 


[aoos  (Ksx  -  vt)  +  ban  (Ksx  -  vt)] 


/ 


The  planetary  and  synoptic  zonal  wind  fields  are: 

ti        1  ^L      A^w      (v     \ 
U  =-- =-— cos  Klx 

f  ay      f  ay     l 

us=4^=4^[a:x)s(Ksx-vt)+bsin(Ksx-vt)j 

f  dy        f  ay 

The  planetary  and  synoptic  meridional  wind  fields  are: 

VL=l_|L=_AWK^sin(KLX) 
I    dx         f 

Vs  =  i4^  =  7NKs[-asin(Ksx-vt)  +  bcos(Ksx-vt)] 
f  dx      f 


(4.21) 


(4.22) 


(4.23) 


(4.24) 


(4.25) 

When  Eqs.  4.17,  4.22  and  4.24  are  substituted  into  Eq.  4.19,  and 
simplified,  the  result  is: 


Q  =  -|rT  P^NflKscosfKLxMKsx-vtK 

r  t  L  oy    ap 

— -N— -Kscos  (KLx)cos  (Ksx  -  vt)+ 

a>'    a? 

W^-KL^sin  (KLx)cos  (Ksx  -  vt)  + 
dy      dP 

W^KL^sin  (KLx)sin  (Ksx  -  vt)  J 


(4.26) 
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Using  the  trigonometric  identities: 

sin(a-p)  =  sinacosp-cosasinp  (4^7) 

cos(a-p)=cosacosp  +  sinasinp  (4  28) 

Eq.  4.26  becomes: 

Q=-^4r-^N^cos(KLx>(^(KSx)cos(vt)-cos(Ksx)sin(vt))+ 
R  f  L   ay     oP 

— — N— Kscos  (KLx)(cos  (Ksx)cos  (vt)  +  sin  (Ksx)sin  (vt)) + 
dy     dP 

W^-KL— sin  (KLx)(cos  (Ksx)cos  (vt)  +  sin  (Ksx)sin  (vt))  + 
ay      dP 

W^KL^sin  (KLx)(sin  (Ksx)cos  (vt)  -  cos  (Ksx)sin  (vt))l 
ay      dP  J 

Using  the  trigonometric  identities: 

sin  a  cos  p  =  i|sin  (a +p )  +  sin  (a  -  B )) 

1  K      H/        y     P//  (4.30) 

cosacosB=l(cos(a+B)+cos(a-B)) 

T  '  "  (4.31) 

sinasinp  =  J4cos  a-p  -cos  a+p)) 

2  (4.32) 

Eq.  4.29  becomes: 


(4.29) 
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Q  =  -  1£A  T  awN||Ks  [cos  (vtXai «KS  +  KL)x)  +  sin «KS  -  KL)x)) 

2  R  f  dy     dP      L 

an(vtKcos((Ks  +  KL)x)  +  cos((Ks-KL)x))l  + 

^N^Ks  [cos  (vt)(cos((Ks  +  KL)x)+cos((Ks  -  KL)x))  + 
dy    dP     L 

sin  (vtKsin  ((Ks  +  KL)x)  +  sin  ((Ks  -  KL)x)) J  + 
W^KL^  [cos  (vtXsin  ((Kl  +  Ks)x)  -  sin  ((Ks  -  KL)x))  + 
an(vtXcos((Ks-KL)x)-cos((Ks  +  KL)x))l  + 
W^KL^  fcos  (vtKcos ((Ks  -  KL)x)  -cos ((Ks  +  KL)x))- 


sin  (vtXsin  ((Kl  +  Ks)x)  -  sin  ((Ks  -  KL)x)) 


]] 


(4.33) 

Eq.  4.33  describes  the  nonlinear  interaction  between  the  planetary-scale  wave 
(wavenumber  four)  and  the  synoptic-scale  wave  (wavenumber  seven).  The 
resulting  waves  are  wavenumber  three  (Ks  -  Kl)  and  wavenumber  11  (Ks  +  Kl). 
Since  synoptic-scale  waves  do  not  propagate  into  the  tropics  (Harris,  1985), 
wavenumber  11  is  ignored  in  Eq.  4.33.  The  resulting  equation,  (4.34),  is  the 
thermal  forcing  function,  Q: 
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Q  =  "  2  R  f  I"  INI*5  [C°S  (vt)sin((Ks  "  K^)-sin(vt)cos((Ks  -  KL)x)J  + 
f^Nf^Ks  ["cos  ( vt)oos  ((Ks  -  KL)x)  +  sin  (vt)sin  ((Ks  -  KLM~|  + 
W^-Kl^  [  cos  (vt)siii ((Ks  -  KL)x)  +  sin  (vt)cos ((Ks  -  KL)x)l  + 

w|Vj^  fcos  (vt)oos  ((KS  -  KL)x)  +  sin  (vt)sin  ((Ks  -  KL)x)l    I 


(4.34) 


2.    Vorticity   Forcing 

The  vorticity  forcing  function  represents  the  nonlinear  advection  of 
vorticity  by  both  the  planetary- scale  and  synoptic-scale  waves.  The  nonlinear 
advection  terms  are: 

A  =  -  Vl-VCs  -  Vs-VCl  (435) 

Eq.  4.35  is  written  in  scalar  form  as  follows: 

aCs     a;s     3sl     3£l 

A  =  -UL— -VL— -Us-^-Vs-^ 

dx  dy  ax  3y  (436) 

When  Eqs.  4.20,  4.21,  4.22,  4.23,  4.24,  and  4.25  are  substituted  into 
Eq.  4.36,  the  result  is: 
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[ 


A  =  A  I  ^WKs 


f 2        ay 


/   ° 


n  „2aN| 


cos(Klx)(-  asin  (Ksx  -  vt) +bcos  (Ksx  -vt))+ 


WkJ  —  -  Ki^  Isin  (KLx)(acos  (Ksx  -  vt)  +  ban  (Ksx  -vt)) 

\ay3      V 

-  K^W  Isin  (KLx)(a20s  (Ksx  -  vt)  +  ban  (Ksx -vt))  ■ 


/   2. 


^-KL 

^y    \ay2 

s| K?— bos  (KLx^- asin  (Ksx- vt)  +  bcos  (Ksx -vt)) 

dy3  dy  / 


] 


(4.37) 

When  like  terms  are  grouped  together  and  the  trigonometric  identities  in 
Eqs.  4.27  and  4.28  are  used  to  simplify  Eq.  4.37,  A  becomes: 
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A  =  A 
f 


[[ 


aw 
ay 


KS 


/   2 


-  KsN 

\ay2  / 


(a3w    ^aw^ 

-NKS K^— 

\  9y3         By  J 


1 


I  cos  (Klx)(-  asin  (Ksxjcos  (vt) +acos  (Ksx)sin  (vt) 
bcos(Ksx)cos(vt)  +  bsin(Ksx)sin(vt))      + 


+ 


I  WKL 


/   3 


3y3       ay    ay 


a*w 


ay: 


KLW 


I] 


sin  (K.Lx)(acos  (Ksx)cos  (vt)  +  asin  (Ksx)sin  (vt)  + 


bsin  (Ksx)cos  (vt)  -  bcos  (Ksx)sin  (vt)") 

-■  (4.38) 

Using  the  trigonometric  identities  of  Eqs.  4.30,  4.31  and  4.32,  Eq.  4.38 


becomes: 
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A=1A 

2f- 


[[ 
[ 


aw 
ay 


Ks 


U2 

'a  n 


\dy2 


-KSN 


a3w    ^aw 


/ 


ay3       dy 


] 


(-  asin  ((Ks  +  KL)x)  - asin  ((Ks  -  KL)x))cos (vt)+ 


(acos((Ks  +  KL)x)  +  acos((Ks  -  KL)x))sin(vt)+ 
(bcos  ((KS  +  KL)x) +bcos  ((Ks  -  KL)x))cos  (vt)+ 


■] 


(bsm((Ks  +  KL)x)  +  bsin((Ks-KL)x))sin(vt)  |   + 


dy 


KL 


\3y2 


KfW 


1 


(asin  ((Ks  +  KL)x)  -  asin  ((Ks  -  KL)x))oos  (vt)  + 


(acos  ((Ks  -  KL)x)  -  acos  ((Ks  +  KL)x))sin  (vt)  + 
(bcos((Ks  -  KL)x)-bcos((Ks  +  KL)x))cos(vt)- 


(bsin  ((KS  +  KL)x)  -  ban  ((Ks  -  KL)x))sin  (vt 


]] 


(4.39) 


Eq.  4.39  describes  the  nonlinear  interaction  between  the  planetary  wave 
(wavenumber  four),  and  the  synoptic  wave  (wavenumber  seven).  The  resulting 
shortwave,  wavenumber  11,  is  ignored.  The  final  form  of  the  vorticity  forcing 
function  is: 
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Maw 
dy 


[[ 
[ 


KS 


^2N    -2J 


w: 


KcN 


NKs 


3  w    2aw 


\9y3         ^ 


I] 


-  asin  ((Ks  -  Kjx)cos  (vt)  +  acos  ((Ks  -  Kjxjsin  (vt)  + 


bcos  ((Ks  -  Kl)x)cos  (vt)  +  ban  ((Ks  -  KL)x)sin  (vt)    I  + 


] 


[ 
[ 


-  rv< 


k3y3         dyj    dy 


ds 


kl 


aV 


By2 


-KLW 


I] 


asin  ((Ks  -  Kjx)oos  (vt)  +  acos  ((Ks  -  KL)x)sin  (vt) 


bcos  ((Ks  -  Kl)x)cos  (vt)  +  bsin  ((Ks  -  Kjx)sin  (vt) 


]] 


(4.40) 
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V.    DESCRIPTION  OF  GENERIC  CYCLONE  MODEL 

The  synoptic-scale  geopotential  height  field  is  given  by  Eq.  4.13,  where  the 
vertical  amplitudes,  a(P,t)  and  b(P,t),  are  obtained  from  the  growth  of  a  generic 
cyclone.  The  intensity  of  the  generic  cyclone  follows  the  development  of 
midlatitude  baroclinic  disturbances  in  the  atmosphere,  using  an  Eady  (1949)  model. 
Eady  assumed  the  atmosphere  was  a  frictionless,  continuously  stratified  fluid 
whose  motion  is  adiabatic,  hydrostatic  and  quasi-geostrophic.  Eady's  original 
study  neglected  the  nonlinear  wave  interactions. 

Peng  (1982)  applied  spectral  methods  to  the  Eady  (1949)  model,  but  included 
frictional  dissipation  and  nonlinear  wave  interactions.  Peng  found  the  behavior  of 
the  baroclinic  wave  depended  on  the  stratification,  S,  frictional  dissipation,  y, 
supercriticality  of  the  vertical  wind  shear,  A,  and  the  fundamental  zonal  wave 
number,  k.  The  Peng  (1982)  model  of  a  baroclinic  wave  is  used  in  this  study  to 
generate  the  generic  cyclone. 

The  atmosphere  is  considered  to  be  an  infinite  channeled  Boussinesq  fluid  with 
a  constant  Brunt-Vaisala  frequency,  Ns.  The  top  and  bottom  boundaries  are  rigid 
walls,  separated  by  a  distance,  D.  The  lateral  boundaries  are  also  rigid  walls, 
separated  by  a  distance,  L.  The  basic  state  is  a  zonal  flow,  U,  with  a  constant 
vertical  wind  shear,  X.  The  quasi-geostrophic  potential  vorticity  equation  that 
describes  the  motion  is: 

(|+Xz|-|q+J(¥,q)=0 

lat     ax/  (5<1) 

where: 
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J(f.g)= *-&-*£ 


dx  dy    dy  3x  (5.2) 


and  the  potential  vorticity  (q)  is: 

q 


i  =  — L+ — L  +  i. — L 


3x2     dy2     S  3Z2  ^ 

with  \\f  being  the  disturbance  streamf unction.  The  vertical  variable  (z)  ranges  from 
minus  1/2  at  the  bottom  of  the  model,  to  1/2  at  the  top. 

The  basic  state  stratification  parameter,  S,  of  the  fluid  is  constant,and  given  by: 

S=' 


NiD2 


*W  (5.4) 

where: 

Ns        -  Brunt-Vaisala  frequency 

D         -  vertical  depth  of  the  fluid 

fo         -  Coriolis  parameter 

L         -  width  of  the  channel 

The  frictional  dissipation,  y,  is  given  by: 

i 

Ro  (5.5) 

where  Ev  is  the  Ekman  number,  and  R0  is  the  Rossby  number. 

The  supercriticality  of  flow,  A,  measures  the  vertical  wind  shear  in  excess  of 
the  critical  vertical  wind  shear  required  for  growth  of  the  baroclinic  wave,  and  is: 
A  =  X-Xc  (5.6) 

where: 
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Ac  = 


4y  (i2tanh2  (I 
k  [(|i  -  tanh  |i )  ( 1  -  juanh  fi)J2 


(5.7) 


represents  the  marginal  shear  required  for  instability  by  linear  theory,  and: 


(5.8) 


where  k  is  the  fundamental  zonal  wave  number.  The  marginal  stability  curve,  as  a 
function  of  (i,  for  S  =  .0628  is  given  in  Fig.  5.1. 


CRITICRL  SHERR 

0.0     2.5     5.0     7.5     10.0 

I    J 

l               i               i 
0.00    0.25       0.50       0.75        1.00        1.25 

MU 

Figure  5.1.    Marginal  Stability  Curve  of  Eady  Model,  Including 

Friction 

The  disturbance  potential  vorticity  (q)  is  equal  to  zero  throughout  the  lifetime 
of  the  cyclone  in  this  model.  Eq.  5.1  simplifies  to: 
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ax2   ay2  s5z2  (59) 

The  nondimensional  boundary  conditions  for  the  model  are: 
-¥=0     aty  =  0,l 

ax 


-.2- 

— -=0    aty  =  0,l 
dydt 


— +  ?iz — — —A.— L+J  \\\f, —  ±Sy 
idt         dx/az       dx       \     3z  / 


(5.10) 


(5.11) 


/   2  2     \ 

a  \i/  a  v 
— -+ — 


=0    atz  =  ±l/2 


W*2     ^i  (5.12) 


where  V  is  the  zonal  average  streamfunction. 

Eqs  5.9,  5.10,  5.11  and  5.12  form  a  closed  system  of  equations.    A  Fourier 

series,  Eq.  5.13,  whose  y  structure  satisfies  Eqs.  5.10  and  5.1 1,  is  used  to  represent 

the  horizontal  streamfunction  field: 

V|/=  X  X  (Cmne"^  +  C^e-"^)  sin  mry  +  £  Q^os  n7^ 

m      n  n  (5.13) 

where  m  and  n  are  the  number  of  zonal  and  meridional  modes,  respectively.   The 
complex  amplitude  coefficient  (Cmn)  of  mode  (m,n),  and  the  complex  conjugate 
(C*mn)  of  Cmn,  are  both  functions  of  time.  Since  only  one  baroclinic  wave  is  used 
in  this  study,  Eq.  5.13  simplifies  to: 
\|/=(Cneikx  +  Cne-ikx)sin7cy  +  Q)icos7r>'  (5.14) 

The  amplitude  coefficients,  Cn,  C*n,  and  C0i  are  determined  by  substituting 
Eq  5.14  into  Eqs.  5.9,  5.10,  5.11  and  5.12.  The  system  of  equations  is  integrated  in 
time  using  a  fourth  order  Runge-Kutta  method  (Peng,  1982).     The  vertical 
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amplitude  coefficients,  a(P,t)  and  b(P,t),  in  Eq.  4.13  are  related  to  the  amplitude 

coefficients  of  Peng  by: 

Cii=acosh(jiz)  +  bsinh(fiz)  (5  15) 

The  variables  a  and  b  are  determined  by  solving  Eq.  5.15  at  z  =  ±  1/2. 

The  nondimensional  coefficients  a  and  b  are  made  dimensional  by  multiplying 
the  them  by  4  000  000  m  (the  characteristic  length  scale  of  the  Eady  model),  40 
m-s-1  (the  characteristic  velocity  scale),  and  .000  1  s_1  (the  Coriolis  parameter  at 

45°  N). 
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VI.    ANALYSIS  OF  RESULTS 

A.   NORMAL  MODE  ANALYSIS 

The  solutions  of  all  three  phases  of  this  study  will  be  analyzed  using  the  vertical 
normal  modes  of  the  primitive  equation  model.  Gill  (1982)  provides  a  good 
general  discussion  of  normal  mode  analysis,  and  Lim  and  Chang  (1987)  derive  the 
vertical  modes  of  a  shallow  water  equation  model  solved  on  an  equatorial  beta- 
plane. 

The  vertical  normal  modes  of  the  primitive  equation  model  used  in  this  study 
are  displayed  in  Figs.  6.1  to  6.8.  There  are  eight  vertical  modes,  consistent  with  the 
eight  levels  of  the  model.  The  first  mode  does  not  have  a  zero-crossing,  and  is 
considered  a  barotropic  mode.  The  second  mode  has  one  zero-crossing.  Since  the 
amplitude  of  the  second  mode  changes  sign  in  the  vertical,  it  is  considered  the 
baroclinic  mode  (Lim  and  Chang,  1987).  The  higher  modes  have  successively 
more  zero-crossings.  The  solutions  to  the  primitive  equation  model  at  each  level 
are  projected  onto  these  vertical  profiles  and  summed  to  produce  the  modal  output. 
The  equivalent  depth  (H)  of  each  vertical  mode  is  provided  in  Table  6.1. 

TABLE  6.1.    EQUIVALENT  DEPTH  OF  VERTICAL  MODES 


Vertical  Mode 

Equivalent  Depth,  m 

1 

7083.1 

2 

209.1 

3 

59.7 

4 

28.4 

5 

16.1 

6 

9.8 

7 

6.2 

8 

3.8 

43 


o.» 


0.7 


<       0$ 

O 

w     o.s 
I 

q 

"^     o^ 


0.3 


0.1 


1 1 1  1 1 1 1 

0  335    0.340        0.345         0.350        0.355        0.360        0.365        0.370         0.375 

AMPLITUDE 


Figure  6.1.     Vertical  Profile  of  Normal  Mode  1 
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Figure  6.2.    Vertical  Profile  of  Normal  Mode  2 
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Figure  6.3.    Vertical  Profile  of  Normal  Mode  3 


Figure  6.4.    Vertical  Profile  of  Normal  Mode  4 
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Figure  6.5.     Vertical  Profile  of  Normal  Mode  5 


Figure  6.6.    Vertical  Profile  of  Normal  Mode  6 
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Figure  6.7.     Vertical  Profile  of  Normal  Mode  7 


Figure  6.8.    Vertical  Profile  of  Normal  Mode  8 
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B.    PHASE  I  RESULTS 

The  forcing  functions  of  Phase  I  will  use  the  development  of  forced  Rossby 
waves  from  Chapter  II.  In  each  experiment  the  forcing  is  applied  between  24°  N 
and  36°  N,  and  the  model  is  integrated  forward  in  time  for  240  hours.  The 
frequency  (co)  and  corresponding  period  of  the  temporal  component  of  the  forcing 
function  (Eq.  4.6)  are  given  in  Table  6.2.  Recall  from  Chapter  II,  that  the  sign  of 
the  variable  a2  determines  whether  the  solution  to  each  vertical  mode  will 
propagate  away  from  the  source  region.  The  variable  a2,  as  a  function  of 
equivalent  depth  and  frequency,  is  given  in  Tables  6.3a  and  6.3b.  The  a2  values  are 
determined  by  first  setting  a  =  K  =  3.797  x  10~7  irr1.  This  value  of  K 
corresponds  to  wavenumber  three  at  30°  N.  The  frequency  is  calculated  for  each 
equivalent  depth  by  using  this  value  of  K  in  Eq.  2.12.  The  variable  a2  is  then 
calculated  for  each  frequency  and  equivalent  depth  using  Eq.  2.1 1 .  Both  fo  and  po 
of  Eqs.  2.1 1  and  2.12  are  calculated  at  30°  N.  When  a2  >  0,  the  solution  for  that 
specific  equivalent  depth  and  frequency  will  propagate  energy  away  from  the 
source  region.  According  to  theory,  only  the  barotropic  mode  (H  =  7083.1  m)  will 
propagate  energy  away  from  the  source  when  the  frequency  (co)  equals 
-206.30  x  10~7  s_1.  As  the  frequency  (co)  decreases,  the  higher  modes  have 
propagating  solutions. 
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TABLE  6.2.    ROSSBY  WAVE  FREQUENCY  AND  PERIOD  OF 

VERTICAL  MODES 


Mode 

Frequency    (co)  ,    s_1 

x   10-7 

Period,  days 

1 

-206.30 

3.5 

2 

-26.13 

27.8 

3 

-8.04 

90.4 

4 

-3.80 

187.4 

5 

-2.22 

327.4 

6 

-1.36 

535.1 

7 

-0.85 

853.2 

8 

-0.53 

1382.0 

TABLE  6.3A.    ALPHA  SQUARED  (M-i  X  10-12)  FOR  HIGH 

FREQUENCY  MODES 


Frequency  (co),  s_1  x  10 

-7 

-206.30        -26.13      -8.04 

-3.88 

7083.1 

0.14            2.60       9.14 

19.18 

209.1 

-2.37            0.14       6.62 

16.66 

Equivalent 

59.7 

-8.85          -6.34       0.14 

10.18 

Depth 

28.4 

-18.89       -16.38     -9.90 

0.14 

(H), 

16.1 

-33.38        -30.86  -24.38 

-14.34 

m 

9.8 

-54.87        -52.36   -45.88 

-35.84 

6.1 

-87.81        -85.29    -78.81 

-68.77 

3.8 

-142.50      -140.00-133.50 

-123.50 

TABLE  6.3B.    ALPHA  SQUARED  (M-i  X  10-12)  FOR  LOW 

FREQUENCY  MODES 


Frequency  (co), 

s-ixl0- 

7 

-2.22 

-1.36        -0.85 

-0.53 

*206.30 

7083.1 

33.67 

55.16        88.10 

142.80 

-0.59 

209.1 

31.15 

52.65        85.58 

140.30 

-3.10 

Equivalent 

59.7 

24.67 

46.17        79.10 

133.80 

-9.58 

Depth 

28.4 

14.63 

36.13        69.06 

123.80 

-19.62 

(H), 

16.1 

0.14 

21.64        54.57 

109.30 

-34.11 

m 

9.8 

-21.35 

0.14        33.08 

87.81 

-55.60 

6.1 

-54.29 

-32.79          0.14 

54.87 

-88.54 

3.8 

-109.00 

-87.52      -54.59 

0.14 

-143.30 

changed  the  sign  of  co  associated  with  mode  1 
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Phase  I  consists  of  seven  experiments.  Each  experiment  will  help  determine 
how  the  primitive  equation  model  reacts  to  the  analytical  forcings  derived  in 
Chapter  II.  The  seven  experiments  are  summarized  in  Table  6.4.  The  first  three 
experiments  will  determine  the  effects  of  a  thermal  forcing  with  different 
equivalent  heights  (H)  or  frequencies  (co).  Experiments  4,  5  and  6  will  determine 
the  effects  of  a  vorticity  forcing  with  different  frequencies.  Since  the  sign  of  co  in 
Experiment  5  is  reversed,  a2  is  negative,  and  all  Rossby  waves  should  be  trapped 
near  the  source.  The  variable  a  in  Experiment  7  is  a  singular  point  in  Eq.  2.15. 
TABLE  6.4.    PHASE  I  EXPERIMENTS 


Experiment 

Thermal  Forcing  Vorticity  Forcing  Equivalent  Depth 

1 

X 

7083.1  m 

2 

X 

209.1  m 

3 

X 

59.7  m 

4 

X 

7083.1  m 

5 

X 

*7083.1  m 

6 

X 

209.1  m 

7 

X 

**209.1  m 

*  sign  of  co  changed 

**a  =  7t/W 

1 .    Experiment  1 

In  this  experiment  the  forcing  function  is  only  added  to  the  thermal 
equation  (Eq.  3.4).  The  mean  or  equivalent  height  (H)  in  Eq.  2.12  is  7083.1  m, 
which  is  exactly  equal  to  the  scale  height  of  the  first  vertical  mode  (the  barotropic 
mode).  The  wave  frequency  (co)  is  -206.30  x  10~7  s_1  and  the  period  is  3.5  days. 
The  wind  velocity  vectors  for  the  lowest  level  (a=  .938)  are  displayed  in  Figs.  6.9 
to  6.18.  The  scale  of  the  wind  vectors  is  shown  in  meters  per  second  at  the  bottom 
of  each  figure. 
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Figure  6.9.    Wind  Vectors  at  24  Hours  for  Phase  I,  Experiment  1 
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Figure  6.10.    Wind  Vectors  at  48  Hours  for  Phase  I,  Experiment  1 
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Figure  6.11.    Wind  Vectors  at  72  Hours  for  Phase  I,  Experiment  1 
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Figure  6.12.    Wind  Vectors  at  96  Hours  for  Phase  I,  Experiment  1 
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Figure  6.13.    Wind  Vectors  at  120  Hours  for  Phase  I,  Experiment  1 
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Figure  6.14.    Wind  Vectors  at  144  Hours  for  Phase  I,  Experiment  1 
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Figure  6.15.    Wind  Vectors  at  168  Hours  for  Phase  I,  Experiment  1 
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Figure  6.16.    Wind  Vectors  at  192  Hours  for  Phase  I,  Experiment  1 
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Figure  6.17.    Wind  Vectors  at  216  Hours  for  Phase  I,  Experiment  1 
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Figure  6.18.    Wind  Vectors  at  240  Hours  for  Phase  I,  Experiment  1 


55 


In  the  first  24  hours  it  is  evident  that  wavenumber  three  develops,  since 
there  is  one  wave  in  each  120°  sector.  The  circulation  is  confined  to  the  latitudes  of 
the  forcing.  By  48  hours  the  overall  magnitude  of  the  vectors  has  increased,  but 
there  still  is  not  significant  propagation  into  the  tropics.  At  the  72  hour  point  the 
magnitude  of  the  wind  vectors  has  decreased,  and  the  circulation  has  pushed 
southward  toward  the  equator.  The  wave  patterns  move  toward  the  west,  but  the 
expected  northeast-southwest  tilt  of  Rossby  waves  is  not  present.  By  96  hours  the 
overall  magnitude  of  the  wind  vectors  is  at  the  minimum,  and  there  is  weak 
circulation  extending  into  the  Southern  Hemisphere.  After  120  hours  the  overall 
magnitude  of  the  wind  vectors  has  increased  to  the  value  at  24  hours,  and  the  weak 
circulation  in  the  tropics  is  not  evident. 

A  similar  pattern  of  changes  occur  at  144,  168,  192,  216  and  240  hours. 
The  three-and-a-half  day  oscillation,  or  beat,  of  the  overall  velocity  vector 
magnitude  is  consistent  with  the  three-and-a-half  day  period  corresponding  to  the 
frequency  of  the  waves  in  this  experiment;  however,  the  beat  frequency  is  not 
consistent  with  the  beta-plane  theory.  The  oscillation  in  magnitude  indicates  the 
wave  pattern  has  not  stabilized,  and  the  time  variation  is  not  purely  sinusoidal,  as 
assumed  in  the  theory . 

The  circulation  patterns  in  Figs.  6.9  to  6.18  do  not  propagate  very  far 
north.  The  trapping  of  waves  to  the  north  may  be  due  to  an  increase  in  the  Coriolis 
parameter  on  the  beta-plane.  As  the  Coriolis  parameter  (f)  increases  in  Eq.  2.1 1, 
a2  becomes  more  negative;  consequently,  the  waves  are  trapped  near  the  source 
latitudes. 

The  wind  vectors  for  the  first  vertical  mode  at  24  and  240  hours  are 
displayed  in  Figs.  6.19  and  6.20.    At  the  24-hour  point,  there  is  significant 
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circulation  in  the  tropics,  but  the  magnitude  is  only  one  meter  per  second  (one  fifth 
of  the  complete  field  shown  in  Fig.  6.9).  This  observation  is  consistent  with  theory, 
in  that  it  is  difficult  for  a  heating  or  mass  source  to  force  the  barotropic  mode  (Lim 
and  Chang,  1987).  After  240  hours,  the  overall  magnitude  of  the  wind  vectors  has 
doubled.  There  is  significant  circulation  south  of  the  equator,  but  the  desired 
northeast-southwest  tilt  of  propagating  Rossby  waves  is  not  evident. 
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Figure  6.19.    Wind  Vectors  at  24  Hours  for  Phase  I,  Experiment  1, 

Mode  1 
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Figure  6.20.    Wind  Vectors  at  240  Hours  for  Phase  I,  Experiment  1, 

Mode  1 

The  pattern  of  the  first  baroclinic  mode  at  24  hours  (Fig.  6.21)  is  similar  to 

the  barotropic  mode  at  24  hours,  but  the  overall  magnitude  of  the  wind  vectors  is 

larger.  This  result  is  similar  to  Lim  and  Chang  (1987).  Lim  and  Chang  found  that 

a  heating  source  is  more  effective  at  forcing  the  higher  order  modes.    The 

circulation  pattern  and  overall  magnitude  of  the  wind  vectors  of  the  first  baroclinic 

mode,  after  240  hours  (Fig.  6.22),  is  similar  to  the  barotropic  mode  at  240  hours, 

except  the  modes  are  1 80°  out  of  phase.  The  overall  magnitude  of  the  wind  vectors 

of  the  higher  modes  (not  shown)  are  less  than  or  equal  to  the  magnitude  of  the 

barotropic  mode  at  24  hours.  According  to  theory,  the  barotropic  mode  is  the  only 

mode  with  a  propagating  solution  (when  H  =  7083.1  m).    The  results  of  this 

experiment  indicate  all  modes  can  propagate  to  some  extent. 
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Figure  6.21.    Wind  Vectors  at  24  Hours  for  Phase  I,  Experiment  1, 

Mode  2 
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Figure  6.22.    Wind  Vectors  at  240  Hours  for  Phase  I,  Experiment  1, 

Mode  2 
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2.    Experiment  2 

In  this  experiment  the  forcing  function  is  again  added  to  the  thermal 
equation,  (Eq.  3.4)  but  the  equivalent  depth  (H)  is  209.09  m.  This  value  of  H 
corresponds  to  the  first  baroclinic  mode,  which  is  mode  two.  The  wave  frequency 
(co)  is  26.13  x  10~7  s_1  and  the  period  is  27.8  days.  The  wind  vector  fields  for  the 
lowest  level  (o  =  .938)  at  24  and  240  hours  are  displayed  in  Figs.  6.23  and  6.24. 
The  patterns  are  not  significantly  different  than  Experiment  1.  The  overall 
amplitude  of  the  wind  vectors  is  much  larger  than  for  Experiment  1.  The  beat 
oscillation  of  the  overall  vector  magnitude  is  not  present,  since  the  period  is  27.8 
days.  The  westward  propagation  is  slower  than  in  Experiment  1,  also  due  to  the 
longer  period.  The  overall  magnitude  of  the  wind  vectors  of  the  barotropic  mode 
at  24  hours  (not  shown)  is  one  meter  per  second,  whereas  the  overall  magnitude  of 
the  baroclinic  wind  vectors  (not  shown)  is  ten  meters  per  second.  The  magnitudes 
of  the  wind  vectors  of  the  higher  modes  are  also  on  the  order  of  the  barotropic 
mode  at  24  hours.  All  the  modes  at  later  times  show  similar  patterns  and 
magnitudes. 


60 


90  N 


60  N 


30°  fT^-A-V 


30"  S 


60°  s  ,- 


90"  S 


j^-grrrrj  ■ 


:  1 


:  ; 


:  :  :  a  :  :  : 


t  t=  i=s: 


:  : 


90  N 


60  N 


^-30  N 


30°  S 


60°S 


120°  W  90°  W    60°  W    30°  W    0° 


-  5. 


30"  E    60°  E    90°  E    120°  E 


90"  S 


Figure  6.23.    Wind  Vectors  at  24  Hours  for  Phase  I,  Experiment  2 
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Figure  6.24.    Wind  Vectors  at  240  Hours  for  Phase  I,  Experiment  2 


61 


3.  Experiment  3 

The  forcing  function  of  this  experiment  is  added  to  the  thermal  equation 
(Eq.  3.4),  and  uses  an  equivalent  depth  (H)  of  59.748  m.  The  frequency  (co)  is 
-8.04  x  10~7  S"1  and  the  period  is  90.4  days.  The  results  are  similar  to  Experiments 
1  and  2,  except  the  overall  magnitude  of  the  velocity  vectors  is  greater  and  the 
westward  propagation  of  the  circulation  patterns  is  slower.  The  magnitude  of  the 
wind  vectors  of  the  barotropic  mode  is  again  much  less  than  the  baroclinic  mode. 

4.  Experiment  4 

In  this  experiment  the  forcing  function  is  only  added  to  the  vorticity 
equation  (Eq.  3.1).  The  equivalent  height  (H)  in  Eq.  2.12  is  7083.1  m,  which  is 
the  height  of  the  barotropic  mode.  The  frequency  (co)  is  -206.30  x  10~7  s_1  and  the 
period  is  3.5  days.  The  wind  vector  fields  for  the  lowest  level  (g  =  .938)  are 
displayed  in  Figs.  6.25  to  6.34.  The  circulation  patterns  are  substantially  different 
than  when  the  source  is  added  to  the  thermal  equations  in  Experiments  1  through  3. 
After  24  hours  there  is  significant  circulation  moving  toward  the  tropics.  By  48 
hours  the  flow  has  crossed  the  equator  and  has  the  northeast-southwest  tilt 
consistent  with  propagating  Rossby  waves.  The  72  through  240  hour  solutions 
show  the  circulation  continues  to  move  toward  the  south.  Notice  again,  that  the 
circulation  does  not  move  too  far  north.  The  overall  magnitude  of  the  vectors  is 
unrealistic  at  the  end  of  the  time  period.  The  large  magnitude  is  caused  by  the 
arbitrary  amplitude  factor  used  in  the  forcing  function.  The  results  of  this 
experiment  are  more  consistent  with  the  midlatitude  beta-plane  theory  than  the 
results  of  Experiments  1  through  3. 
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Figure  6.25.    Wind  Vectors  at  24  Hours  for  Phase  I,  Experiment  4 
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Figure  6.26.    Wind  Vectors  at  48  Hours  for  Phase  I,  Experiment  4 
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Figure  6.27.    Wind  Vectors  at  72  Hours  for  Phase  I,  Experiment  4 


90    N 


60°  N 


30°  N 


90°  S 


90"  N 


120    W      90    W 


60    W 


30UE 

20.0 


90   E 


120    E 


Figure  6.28.    Wind  Vectors  at  96  Hours  for  Phase  I,  Experiment  4 


64 


90  N 


! 


i  - 


im 


111' 


^  B  &  4- 


Z'-A 


r-  -  <  ;  V^^^^-A 

F,te 


=  22 

§  §  =  ^ 


:  H  M 


II M 


Itti 


■*## 


FT 


120°  W  90°  W    60°  W    30°  W 


90UN 


60  N 


30"  N 


-  20.0 


60"  E 


90  E    120  E 


Figure  6.29.    Wind  Vectors  at  120  Hours  for  Phase  I,  Experiment  4 
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Figure  6.30.    Wind  Vectors  at  144  Hours  for  Phase  I,  Experiment  4 
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Figure  6.31.    Wind  Vectors  at  168  Hours  for  Phase  I,  Experiment  4 
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Figure  6.32.    Wind  Vectors  at  192  Hours  for  Phase  I,  Experiment  4 
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Figure  6.33.    Wind  Vectors  at  216  Hours  for  Phase  I,  Experiment  4 
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Figure  6.34.    Wind  Vectors  at  240  Hours  for  Phase  I,  Experiment  4 
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The  overall  magnitude  of  the  wind  vectors  of  the  barotropic  mode  at  24 
hours  (Fig.  6.35)  is  five  times  larger  than  the  magnitude  of  the  wind  vectors  of  the 
baroclinic  mode  at  the  same  time  (Fig.  6.36).  The  vorticity  source  is  much  more 
efficient  at  forcing  the  barotropic  mode. 
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Figure  6.35.    Wind  Vectors  at  24  Hours  for  Phase  I,  Experiment  4, 

Mode  1 
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Figure  6.36.    Wind  Vectors  at  24  Hours  for  Phase  I,  Experiment  4, 

Mode  2 

5.    Experiment  5 

This  experiment  is  exactly  the  same  as  Experiment  4,  except  the  sign  of  the 

frequency  (co)  in  Eq.  2.12  is  changed.    With  a  positive  frequency  (w),  \i  >  0  in 

Eq.  2.26  and  the  circulation  should  be  trapped  near  the  source.  The  wind  velocity 

vectors  for  the  lowest  level  (a  =  .938)  at  24  and  240  hours  are  displayed  in 

Figs.  6.37  and  6.38.  The  overall  magnitude  of  the  vectors  is  significantly  smaller 

than  those  of  Experiment  4.    However,  the  circulation  did  propagate  past  the 

equator,  but  to  a  much  less  extent.    The  propagation  is  most  likely  due  to  fast 

moving  gravity  waves.  The  modal  analysis  is  consistent  with  the  complete  solution, 

and  the  results  of  Experiment  4. 
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Figure  6.37.    Wind  Vectors  at  24  Hours  for  Phase  I,  Experiment  5 
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Figure  6.38.    Wind  Vectors  at  240  Hours  for  Phase  I,  Experiment  5 
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6.    Experiment  6 

In  this  experiment  the  forcing  is  only  added  to  the  vorticity  equation 
(Eq.  3.1).  The  equivalent  height  (H)  is  209.09  m,  which  is  the  height  of  the  first 
baroclinic  mode.  The  frequency  (co)  is  26.13  x  10"7  s-1  and  the  period  is  27.8  days. 
The  wind  velocity  vectors  for  the  lowest  level  (o  =  .938)  at  24  and  240  hours  are 
displayed  in  Figs.  6.39  and  6.40.  The  circulation  patterns  show  a  significant  tilt 
toward  the  southwest,  but  the  circulation  does  not  propagate  significantly  into  the 
Southern  Hemisphere.  The  overall  magnitude  of  the  barotropic  mode  vectors  (not 
shown)  is  again  much  larger  than  the  baroclinic  mode  vectors  (not  shown)  at  the 
24-hour  point. 
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Figure  6.39.    Wind  Vectors  at  24  Hours  for  Phase  I,  Experiment  6 
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Figure  6.40.    Wind  Vectors  at  240  Hours  for  Phase  I,  Experiment  6 
7.    Experiment  7 

In  this  experiment  the  forcing  function  is  added  to  the  vorticity  equation 
(Eq.  3.1).  The  equivalent  height  (H)  in  Eq.  2.12  is  209.09  m,  which  is  exactly 
equal  to  the  first  baroclinic  mode.  The  variable  a  equals  n  I W,  where  W  is  half  the 
width  of  the  forcing  function  in  degrees  latitude.  This  value  of  a  is  a  singular  point 
in  the  complete  propagating  solution  (Eqs.  2.17  and  2.18).  The  wind  velocity 
vectors  for  the  lowest  level  (a  =  .938)  are  displayed  in  Figs.  6.41  to  6.50. 
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Figure  6.41.    Wind  Vectors  at  24  Hours  for  Phase  I,  Experiment  7 
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Figure  6.42.    Wind  Vectors  at  48  Hours  for  Phase  I,  Experiment  7 
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Figure  6.43.    Wind  Vectors  at  72  Hours  for  Phase  I,  Experiment  7 
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Figure  6.44.    Wind  Vectors  at  96  Hours  for  Phase  I,  Experiment  7 
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Figure  6.45.    Wind  Vectors  at  120  Hours  for  Phase  I,  Experiment  7 
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Figure  6.46.    Wind  Vectors  at  144  Hours  for  Phase  I,  Experiment  7 
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Figure  6.47.    Wind  Vectors  at  168  Hours  for  Phase  I,  Experiment  7 
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Figure  6.48.    Wind  Vectors  at  192  Hours  for  Phase  I,  Experiment  7 
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Figure  6.49.    Wind  Vectors  at  216  Hours  for  Phase  I,  Experiment  7 
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Figure  6.50.    Wind  Vectors  at  240  Hours  for  Phase  I,  Experiment  7 
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The  overall  magnitude  of  the  wind  vectors  is  approximately  the  same  as  in 
Experiment  6,  but  the  circulation  pattern  has  more  of  a  northeast-southwest  tilt. 
The  frequency  (co)  of  this  forcing  function  is  -3.05  x  lO7  s_1  and  the  period  is 
238.6  days.  The  westward  movement  of  the  circulation  pattern  is  much  slower  than 
in  Experiment  6,  due  to  the  longer  period.  The  modal  analysis  is  consistent  with  the 
analyses  of  Experiments  4  through  6. 
8.    Phase  I  Summary 

The  results  of  Phase  I  indicate  that  forced  Rossby  wave  forcings  applied  to 
the  thermal  equation  do  not  produce  significant  Rossby-type  cold  surge  responses. 
However,  when  the  same  forcings  are  applied  to  the  vorticity  equation,  the 
responses  are  consistent  with  Rossby-beta-plane  theory.  Forcing  the  thermal 
equation  consistently  produces  a  stronger  response  in  the  baroclinic  mode,  and 
forcing  the  vorticity  equation  produces  a  stronger  barotropic  response.  When  the 
sign  of  the  frequency  (co)  is  changed,  the  solution  does  propagate  southward.  The 
southward  propagating  solution  is  not  consistent  with  Rossby  wave  theory.  The 
propagating  waves  in  this  case  could  be  gravity  waves  since  the  overall  magnitude 
of  the  velocity  vectors  is  small. 

C.    PHASE  II  RESULTS 

The  forcing  functions  of  Phase  II  use  the  impulse  function  described  in 
Chapter  V.  The  maximum  amplitude  of  the  forcing  is  at  24  hours  and  a  =  .375.  In 
each  experiment  the  forcing  is  applied  between  24°  N  and  36°  N,  and  the  model  is 
again  integrated  forward  in  time  for  240  hours.  The  maximum  amplitude  of  the 
forcing  occurs  at  24  hours  (as  shown  in  Fig.  4.3).  The  two  experiments  of  Phase  II 
are  summarized  in  Table  6.5.  Experiments  8  and  9  will  compare  the  effects  of  the 
impulse  forcing  in  either  the  thermal  or  vorticity  equations. 


78 


TABLE  6.5. 

PHASE  II  EXPERIMENTS 

Experiment 

Thermal  Forcing  Vorticity  Forcing 

8 
9 

X 

X 

1.    Experiment  8 

In  this  experiment  the  impulse  source  is  only  added  to  the  thermal  equation 
(Eq.  3.4).  The  wind  velocity  vector  fields  for  the  lowest  level  (a  =  .938)  are 
displayed  in  Figs.  6.51  through  6.60.  After  24  hours  a  wavenumber  three 
circulation  develops  near  the  source  latitudes.  The  overall  magnitude  of  the 
velocity  vectors  is  relatively  small  and  the  circulation  does  not  propagate  very  far 
south.  At  48  hours  the  magnitude  increases  and  the  circulation  pattern  moves 
toward  the  west,  but  the  circulation  is  still  confined  to  the  source  latitudes.  By  144 
hours  there  is  a  slight  northeasterly  flow  toward  the  equator.  The  magnitude  of  the 
northeasterly  flow  does  not  significantly  increase  after  240  hours. 
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Figure  6.51.    Wind  Vectors  at  24  Hours  for  Phase  II,  Experiment  8 
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Figure  6.52.    Wind  Vectors  at  48  Hours  for  Phase  II,  Experiment  8 
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Figure  6.53.    Wind  Vectors  at  72  Hours  for  Phase  II,  Experiment  8 
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Figure  6.54.    Wind  Vectors  at  96  Hours  for  Phase  II,  Experiment  8 
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Figure  6.55.    Wind  Vectors  at  120  Hours  for  Phase  II,  Experiment  8 
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Figure  6.56.    Wind  Vectors  at  144  Hours  for  Phase  II,  Experiment  8 
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Figure  6.57.    Wind  Vectors  at  168  Hours  for  Phase  II,  Experiment  8 
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Figure  6.58.    Wind  Vectors  at  192  Hours  for  Phase  II,  Experiment  8 
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Figure  6.59.    Wind  Vectors  at  216  Hours  for  Phase  II,  Experiment  8 
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Figure  6.60.    Wind  Vectors  at  240  Hours  for  Phase  II,  Experiment  8 

The  overall  magnitude  of  the  velocity  vectors  of  the  barotropic  mode  (not 
shown)  is  on  the  order  of  0.5  m-s-1.  The  overall  magnitude  of  the  first  baroclinic 
mode  velocity  vectors  (not  shown)  is  ten  times  greater.  The  larger  baroclinic 
response  is  consistent  with  Phase  I.  The  circulation  patterns  of  the  modal  analyses 
are  not  significantly  different  than  the  total  patterns  shown  in  Figs.  6.51  through 
6.60. 

2.    Experiment  9 

In  this  experiment  the  impulse  function  is  only  added  to  the  vorticity 
equation  (Eq.  3.1).  The  wind  velocity  vector  fields  for  the  lowest  level  (a=  .938) 
are  displayed  in  Figs.  6.61  through  6.70.  After  24  hours  a  wavenumber  three 
circulation  develops  and  significant  circulation  patterns  propagate  toward  the 
equator.  By  48  hours  the  magnitude  doubles  and  the  northeasterly  flow  intensifies. 
By  72  hours  there  is  significant  flow  past  the  equator.    A  definite  northeast- 
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southwest  tilt  in  the  circulation  pattern  is  evident  at  96  hours.  This  tilt  is  much 
more  pronounced  than  in  Experiment  8.  By  192  hours  there  are  closed  circulation 
cells  in  the  tropics.  It  is  interesting  to  note  that  the  circulation  patterns  of  this 
experiment  also  propagate  north  of  the  forcing  region. 
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Figure  6.61.    Wind  Vectors  at  24  Hours  for  Phase  II,  Experiment  9 
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Figure  6.62.    Wind  Vectors  at  48  Hours  for  Phase  II,  Experiment  9 
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Figure  6.63.    Wind  Vectors  at  72  Hours  for  Phase  II,  Experiment  9 
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Figure  6.64.    Wind  Vectors  at  96  Hours  for  Phase  II,  Experiment  9 
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Figure  6.65.    Wind  Vectors  at  120  Hours  for  Phase  II,  Experiment  9 
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Figure  6.66.    Wind  Vectors  at  144  Hours  for  Phase  II,  Experiment  9 
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Figure  6.67.    Wind  Vectors  at  168  Hours  for  Phase  II,  Experiment  9 
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Figure  6.68.    Wind  Vectors  at  192  Hours  for  Phase  II,  Experiment  9 
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Figure  6.69.    Wind  Vectors  at  216  Hours  for  Phase  II,  Experiment  9 
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Figure  6.70.    Wind  Vectors  at  240  Hours  for  Phase  II,  Experiment  9 

The  velocity  vector  patterns  of  the  first  and  second  modes  are  similar  to 
the  complete  vector  patterns.    The  overall  magnitude  of  the  barotropic  mode 
vectors  is  ten  times  greater  than  the  baroclinic  mode  vectors  at  24  hours.    The 
larger  response  of  the  barotropic  mode  is  consistent  with  the  results  of  Phase  I. 
3.    Phase  II  Summary 

The  results  of  Phase  II  indicate  that  the  thermal  impulse  forcing  is  not  an 
efficient  method  of  producing  a  cold  surge  response.  After  240  hours  a  small  but 
significant  northeasterly  flow  develops,  and  circulation  patterns  are  evident  near 
the  equator;  however,  ten  days  is  too  slow  to  be  realistic.  The  impulse  forcing 
applied  to  the  vorticity  equation  does  produce  a  significant,  timely  response.  The 
northeast-southwest  tilt  in  the  circulation  patterns  is  similar  to  Bashford  (1985)  and 
Lim  and  Chang  (1981).  The  results  are  consistent  with  Phase  I,  in  that  the  vorticity 
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forcing  is  more  efficient  in  exciting  a  barotropic  response,  and  the  thermal  forcing 
a  baroclinic  response. 

D.   PHASE  III  RESULTS 

The  forcing  functions  of  Phase  III  use  the  generic  cyclone  described  in 
Chapter  V.  In  each  experiment  the  baroclinic  zone  is  applied  between  30°  N  and 
60°  N.  This  latitudinal  range  is  representative  of  a  synoptic  wave  moving  through 
the  planetary  wave  over  East  Asia.  The  latitudinal  structure  of  both  the  thermal 
and  vorticity  forcing  functions  are  shown  in  Fig.  6.71.  The  magnitude  of  each 
forcing  equals  zero  in  the  Southern  Hemisphere  and  outside  the  source  region.  The 
model  is  integrated  forward  in  time  for  240  hours.  The  maximum  amplitude  of  the 
forcings  occur  at  the  24-hour  and  216-hour  points  (Fig.  6.72).  In  the  vertical,  the 
maximum  amplitude  of  both  thermal  (Fig.  6.73)  and  vorticity  forcings  (Fig.  6.74) 
occur  at  the  surface.  The  three  experiments  in  Phase  III  will  determine  which 
forcing  is  dominant  in  the  generation  of  cold  surges.  The  experiments  are 
summarized  in  Table  6.6. 

TABLE  6.6.  PHASE  III  EXPERIMENTS 


Experiment Thermal  Forcing  Vorticity  Forcing 

10  XX 

11  X 

12  X 
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Figure  6.71.    Latitudinal  Structure  of  Phase  III  Forcings 
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Figure  6.72.    Temporal  Structure  of  Phase  III  Forcings 
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1.    Experiment   10 

In  this  experiment  the  forcing  is  applied  to  both  the  thermal  (Eq.  3.4)  and 
vorticity  (Eq.  3.1)  equations.  The  wind  velocity  vector  fields  for  the  lowest  level 
(a  =  .938)  are  displayed  in  Figs.  6.75  through  6.84.  After  24  hours  a  wavenumber 
three  circulation  develops  near  the  source  latitudes  and  there  is  some  propagation 
toward  the  south.  By  48  hours  a  five  meters  per  second  northeasterly  wind,  typical 
of  an  East  Asian  cold  surge,  is  evident  to  approximately  10°  N.  After  72  hours  the 
northeasterly  flow  extends  past  the  equator  into  the  Southern  Hemisphere.  The 
northeast-southwest  tilt  in  the  circulation  south  of  the  source  latitudes,  and 
northwest-southeast  tilt  north  of  the  source  is  evident  at  this  time.  The  tilted 
circulation  pattern  is  characteristic  of  propagating  Rossby  waves.  The  propagating 
wave  continues  in  this  manner,  and  at  240  hours  there  is  identifiable  circulation  at 
30°  S. 
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Figure  6.75.    Wind  Vectors  at  24  Hours  for  Phase  III,  Experiment  10 
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Figure  6.76.    Wind  Vectors  at  48  Hours  for  Phase  III,  Experiment  10 


90    N 


60UK 


30°  N    s 


30US    »- 


60    S    i 


90"  S 


wm 


,i.i      r    —   <?  s>    >> 


1 


_. r»-_ 


r-    ~»    -»    K     *      i     * 

A  A  i    i. 


-.  -.  ">    > 

•30    N 


90    N 


30"  S 


60"  S 


120°  W     90°  W         60°  W         30°  W  0C 


90°  S 


30    E  60    E  90   E  120    E 

-    20.0 


Figure  6.77.    Wind  Vectors  at  72  Hours  for  Phase  III,  Experiment  10 
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Figure  6.78.    Wind  Vectors  at  96  Hours  for  Phase  III,  Experiment  10 
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Figure  6.79.    Wind  Vectors  at  120  Hours  for  Phase  III,  Experiment  10 
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Figure  6.80.    Wind  Vectors  at  144  Hours  for  Phase  III,  Experiment  10 
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Figure  6.81.    Wind  Vectors  at  168  Hours  for  Phase  III,  Experiment  10 
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Figure  6.82.    Wind  Vectors  at  192  Hours  for  Phase  III,  Experiment  10 
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Figure  6.83.    Wind  Vectors  at  216  Hours  for  Phase  III,  Experiment  10 
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Figure  6.84.    Wind  Vectors  at  240  Hours  for  Phase  III,  Experiment  10 

The  barotropic  (Fig.  6.85)  and  baroclinic  (Fig.  6.86)  mode  analyses  at  24 
hours  show  similar  circulation  patterns  as  the  complete  solution  at  24  hours.  The 
overall  magnitude  of  the  barotropic  mode  vector  field  is  five  times  larger  than  the 
baroclinic  vector  field.  The  larger  barotropic  response  suggests  the  vorticity 
forcing  is  dominant  in  the  generic  cyclone.  The  higher  modes  have  similar 
circulation  patterns,  but  the  overall  magnitude  of  the  vector  field  is  much  smaller. 
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Figure  6.85.    Wind  Vectors  at  24  Hours  for  Phase  III,  Experiment  10, 

Mode  1 


90    N 


60    N 


30    N  j— f 


LfJJJji    v^jJ-A  U£/jj   tv^il 


30*  S 


60"  S 


90°  S 


90"  N 


60  N 


30*  N 


30  S 


60°  S 


120  W  90  \V    60  W    30  W    0 


30°  E    60°  E    90CE    120°  E 


90°  S 


Figure  6.86.    Wind  Vectors  at  24  Hours  for  Phase  III,  Experiment  10, 
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2.    Experiment   11 

In  this  experiment  the  forcing  is  only  applied  to  the  thermal  equation 
(Eq.  3.4).  After  24  hours  a  wavenumber  three  circulation  develops  in  the  lowest 
level  (a  =  .938),  but  the  circulation  is  rather  chaotic  (Fig.  6.87).  The  circulation 
remains  trapped  for  the  whole  time  period.  After  240  hours  (Fig.  6.88)  the 
sinusoidal  circulation  pattern  is  still  within  the  source  latitudes.  It  is  apparent  that 
the  thermal  forcing  does  not  play  a  significant  role  in  producing  cold  surges  in  this 
model.  The  overall  magnitude  of  the  baroclinic  mode  vector  field  at  24  hours  (not 
shown)  is  four  times  larger  than  the  magnitude  of  the  barotropic  field.  The  larger 
baroclinic  response  is  similar  to  the  results  of  Phase  I  and  II. 
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Figure  6.87.    Wind  Vectors  at  24  Hours  for  Phase  III,  Experiment  11 
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Figure  6.88.    Wind  Vectors  at  240  Hours  for  Phase  III,  Experiment  11 

3.  Experiment   12 

In  this  experiment  the  forcing  is  only  applied  to  the  vorticity  equation 
(Eq.  3.1).  Since  the  vorticity  forcing  is  dominant,  and  the  thermal  forcing  does 
not  propagate  away  from  the  source,  the  circulation  patterns  and  magnitudes  are 
not  significantly  different  than  the  combined  thermal  and  vorticity  forcing  of 
Experiment  11. 

4.  Phase  III  Summary 

The  results  of  Phase  III  indicate  the  vorticity  forcing  is  the  dominant 
component  in  the  nonlinear  interaction  of  a  planetary  wave  and  a  synoptic-scale 
generic  cyclone.  The  thermal-only  forcing  does  not  propagate  away  from  the 
source  and  is  less  intense  than  the  vorticity-only  forcing.  The  vorticity  forcing 
(and  thus  the  combined  thermal  and  vorticity  forcing)  does  produce  a  northeasterly 
wind  within  48  hours,  which  is  typical  of  East  Asian  cold  surges.  The  northeasterly 
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flow  continues  into  the  Southern  Hemisphere,  which  suggests  a  tropical  response  to 
a  midlatitude  forcing  (Lim  and  Chang,  1987). 
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VII.    CONCLUSION 

The  objective  of  this  study  was  to  determine  whether  the  nonlinear  interaction 
between  a  planetary-scale  wave  (wavenumber  four)  and  a  rapidly  growing 
synoptic-scale  wave  (wavenumber  seven)  could  force  a  northeasterly  cold  surge. 
The  nonlinear  interaction  of  these  waves  produced  a  wavenumber  three  forcing. 
The  basic  model  was  a  linearized  global  spectral  model  with  eight  levels.  Fourteen 
experiments  with  three  different  types  of  source  functions  were  applied  to  the 
thermal  and/or  vorticity  equations.  The  solutions  of  each  experiment  were 
projected  onto  the  vertical  normal  modes  of  the  primitive  equation  model. 

The  results  of  each  phase  of  the  experiments  indicate  a  midlatitude  thermal 
forcing  does  not  produce  a  significant  cold  surge  response.  The  impulse  forcing, 
similar  to  Lim  and  Chang  (1981),  applied  to  the  thermal  equation  produced  a  weak 
cold  surge  response  after  240  hours;  however,  ten  days  is  longer  than  the  typical 
two-day  surge  generation  period  observed  in  the  East  Asian  winter  monsoon.  The 
difference  between  the  equatorial  beta-plane  model  used  by  Lim  and  Chang  (1981), 
and  the  spherical  primitive  equation  model  used  in  this  study  will  require  further 
explanation.  However,  the  thermal  forcing  consistently  produced  a  stronger 
response  in  the  baroclinic  mode. 

The  midlatitude  vorticity  forcing,  on  the  other  hand,  produced  a  stronger 
response  in  the  barotropic  mode.  The  wavenumber  three  vorticity  forcing 
produced  a  significant  cold  surge  response  within  48  hours.  The  cold  surge 
produced  by  the  vorticity  forcing  is  consistent  with  the  observed  cold  surge  in  the 
East  Asian  winter  monsoon. 
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The  lack  of  cold  surge  response  to  a  midlatitude  thermal  forcing,  and  the 
significant  response  to  a  midlatitude  vorticity  forcing  may  be  explained  by 
geostrophic  adjustment  theory.  In  midlatitude  geostrophic  adjustment  the  wind 
field  tends  to  adjust  to  the  mass  field.  Thus,  when  the  forcing  is  prescribed  by 
heating,  the  local  adjustment  (to  rotational  motion)  is  more  complete.  When  the 
forcing  is  prescribed  by  vorticity,  the  adjustment  will  be  less  efficient  and  a 
considerable  amount  of  energy  will  propagate  away  in  divergent  motion.  The 
equatorward  propagation  studied  by  Lim  and  Chang's  (1981)  theory  is  through 
equatorial  Rossby  waves  that  have  a  significant  divergent  component.  Thus,  the 
divergent  motion  may  offer  a  possible  explanation  for  the  different  responses  to 
the  two  types  of  forcings. 

The  initial  condition  of  this  study  did  not  include  any  basic  current  (mean 
flow).  In  reality,  the  cold  surge  is  superimposed  on  an  existing  northeasterly 
monsoonal  flow  at  the  surface.  It  would  be  interesting  to  see  if  the  presence  of  an 
initial  mean  flow  would  alter  the  results  of  the  forcings.  Also,  the  nonlinear  effects 
of  topographic  channeling  play  an  important  role  in  the  cold  surge  response.  Using 
some  type  of  topographic  barrier  may  lead  to  interesting  conclusions. 
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